Modulation and Measurement of Time-Energy Entangled Photons 
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We describe a proof-of-principal experiment demonstrating a Fourier technique for measuring the 
shape of biphoton wavepackets. The technique is based on the use of synchronously driven fast 
modulators and slow (integrating) detectors. 

PACS numbers: 42.50.Gy, 32.80.Qk, 42.50.Ex, 42.65.Lm 
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The characterization of the joint quantum state of 
time-energy entangled photons is a long-standing prob- 
lem of interest to workers in quantum information 
processing[l[ This interest is. in part, motivated by 
the quantum property where either dispersion Q or mod- 
ulation ^ that is experienced by a photon in one channel 
may be negated by the dispersion or modulation that is 
experienced by a second photon in a different channel. By 
using appropriate spectral tools, the measurement of fre- 
quency domain correlation is readily accomplished. The 
measurement of temporal correlation is more difficult: If 
the temporal resolution of the photon detector is high 
enough, then a time-to-digital converter (TDC) may be 
used to directly measure it 0, @|. When the detectors 
do not have sufficient speed the measurement alterna- 
tives are limited. The technique of Hong-Ou-Mandel 
docs not respond to even order dispersion Q and, though 
often used for this purpose, only measures the bipho- 
ton waveform when the waveform is transform limited. 
When sufficient count rate is available, sum frequency 
generation in a nonlinear crystal may be used as an ultra- 
fast correlator [l^l . If the down-converting crystal is 
synchronously pumped, then femtosecond time-scale up- 
conversion may be used for temporal correlation . 
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FIG. 1: Measurement of biphoton wave functions. Ampli- 
tude modulators in the signal and idler paths are driven syn- 
chronously. The coincidence count rate versus frequency is 
measured by slow detectors. 

This Letter describes a proof-of-principal experiment 
that demonstrates a new approach to the problem of 
measurement of biphoton wavefunctions. Fig. 1 shows 
the essential idea. Spontaneously generated signal and 
idler photons, or as in the experimental portion of this 
work. Stokes and anti-Stokes photons, are incident on 
synchronously driven sinusoidal amplitude modulators. 
The coincidence count rate between single photon count- 



ing modules (SPCMs) is measured as a function of the 
sinusoidal modulation frequency. The SPCMs are slow as 
compared to the pulse width of the biphoton wave func- 
tion, so that they integrate over the wavefunction. With 
T equal to the relative arrival time of the signal and idler 
photons, the inverse Fourier transform of the (frequency 
domain) measurement of coincidence counts versus fre- 
quency yields the Glauber correlation function G^(r) and 
therefore the square of the absolute value of the bipho- 
ton wavefunction By definition, time to frequency 
Fourier transformation is the result of multiplication by 
a sinusoidal function and integration over time. Here, 
this operation is accomplished by modulation (multipli- 
cation) followed by slow detection (integration). 

In recent experiments, Kolchin et al. [l^ have demon- 
strated electro-optic modulation of single photons by us- 
ing one photon of a biphoton pair to set the time origin 
for the modulation of the second photon. Here, modu- 
lation of the biphoton wavefunction is attained without 
establishing a time origin for the modulator; i.e., the ab- 
solute common phase of the sinusoidal modulating signal 
is of no consequence. Entangled signal and idler photons 
may be thought of as arriving at the modulators at the 
same, but random, time. If the modulators are driven 
at different frequencies, or one channel is modulated and 
the other is not, then the Glauber correlation function is 
not changed by the modulators. If the modulators are 
driven at the same frequency, then, as will be shown be- 
low, the coincidence count rate is modulated. We note 
that some aspects of modulation of biphotons have been 
considered by Tsang et al. [l3| . 

Theory: We work in the Heisenberg picture with ex- 
pectations taken against the vacuum state at the input of 
the parametric generator [14|. We let asfi{t) and aifi{t) 
denote the signal and idler operators at the output of the 
parametric generator and before the amplitude modula- 
tors. With t and t + t denoting the arrival times of the 
signal and idler photons at their respective detectors, the 
unmodulated Glauber intensity correlation function is 

then G[)^^(r) = (aJ_o(t)aJ o(i + T)a,,o(t + r)a»,o W)- When 
amplitude modulators mi{t) and m2{t) are inserted be- 
tween the downconversion source and the detectors, the 
field operators are modulated so that aa{t) = ■mi{t)as,o{t) 
and ai(t) = m2{t)ai_Q(t). The modulated correlation 
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function is then 



(1) 



For a source with time independent statistics the un- 

(2) 

modulated correlation function Gq '{t) is a function of 
the difference of the detection times of the two photons r 
and does not depend explicitly on the arrival time t of the 
first photon. But as a result of the random time of arrival 
of the photon pair with regard to the phased modulators, 
the correlation function of the modulated biphoton does 
depend on both t and r. We average G ^[ {t^ t -j- t) over 
a period T of the modulating frequency to form the time 
averaged correlation function 
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Combining Eq. (1) and Eq.(2), we obtain 



M{t) = 



M{t)G^^\t) 
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T)\^dt. 
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where M.{t) is the intensity correlation function of the 
modulators in the signal and idler channels. We thus 
have the result that the time-averaged Glauber correla- 
tion function of the modulated biphoton wavefunction is 
equal to the correlation function of the modulators mul- 
tiplied by the unmodulated Glauber correlation function. 

If one channel is modulated and the other is not, 
Ai{T) = 1/2, and the intensity correlation function is 
not modulated. If both channels are modulated by si- 
nusoidal amplitude modulators with frequency uj and 
a common phase ip, that is by modulators mi(t) = 
7712(0 — coslujt -f- if), then irrespective of this phase, 
M{t) = 1/4+ 1/8 cos(2wr). 

The essence of the Fourier technique is now apparent: 
If the detectors are slow in the sense that they inte- 
grate the coincidence counts for a time that is long as 
compared the the length of the biphoton wave packet 
, but short as compared to the inverse rate of bipho- 
ton generation, then the integrated coincidence count 



is J,G\nr)dr = 1/8 /„ [2 + cos(2c.r)] G{,^^ (r)dr. We 
neglect the DC term, and normalize to obtain the Fourier 
cosine transform pair 
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G^q\t) cos{2uJT)dT 
F{2uj) cos(2ujT)duj (4) 



With the unmodulated biphoton wavefunction denoted 
by 0(t), and taking the rate of generated paired photons 
as R, then as shown by Rubin et al. [15|] and Kolchin 
[1^, G[,^^(r) = |0(r)|2 + R^. Assuming a bin width 



5T, and photodetectors having a quantum efficiency of 
e, the coincidence count rate Rc{t) is related to Gq (t) 



hy R,{t) = {eHT)Gf\T) 
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FIG. 2: (color online) (a) Schematic of experiment. Long 
biphotons are generated in the cold-atom Rb cell. Using an 
optical fiber, the anti-Stokes photon is delayed by 175 ns and 
the photons are modulated by synchronously driven sinusoidal 
modulators (b) Energy level diagram for paired photon gen- 
eration in Rb. 

Experiment: Our experiment makes use of long bipho- 
tons that are produced using the techniques of electro- 
magnetically induced transparency and slow light. The 
use of long biphotons allows us to compare the correla- 
tion function measured by our Fourier transform tech- 
nique with a direct measurement using fast detectors 
and a TDC. The experimental configuration is shown 
in Fig. 2 (a). Paired photons are generated with cold 
Rb atoms using the method of Balic et al. [1]. We apply 
strong counterpropagating pump and coupling lasers (not 
shown) to produce phase matched counter-propagating 
pairs of time-energy entangled Stokes and anti-Stokes 
photons. We use a ^^Rb two-dimensional magneto-optic 
trap with an optical depth which can be varied between 
10 to 60 17| to generate biphotons with temporal lengths 
between 50 and 900 ns. The inset in Fig. 3(a) shows the 
biphoton wavefunction obtained at an optical depth of 
35 as directly measured using a TDC. Two features are 
of interest. First, the width of the wavefunction is de- 
termined by the slow group velocity of the anti-Stokes 
photon and varies linearly with the optical depth 



Second, the distinctive sharp feature at the leading edge 
is a Sommerfield-Brillouin precursor that ensures that 
the earliest signal reaches a detector at the speed of light 
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The generated Stokes and anti-Stokes photons are 
transmitted through 10 GHz electro-optic amplitude 
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FIG. 3: (color online) Modulation of the biphoton correlation 
function, (a) mi{t) = 1 and m2(i) = cos(ti;t + (p). The inset 
in (a) shows the correlation function with both modulators 
open, (b) mi{t) = m^it) = cos{ujt + if). Here uj — 2nx 35 
xlO®. 



FIG. 4: (color online) Fourier Transform measurement tech- 
nique at an optical depth of 15. (a) Frequency domain data, 
(b) Fourier Cosine transform (gray), and temporal correla- 
tion (black). The term modulation frequency on the x-axis of 
(a) refers to the applied frequency. The observed modulation 
frequency [Eq.(4)] is a factor of 2 higher. 



modulators (Eospacc Inc.) with a half-wave voltage, 14^ 
of 1.3V. To obtain a perfectly sinusoidal output, the mod- 
ulators arc biased at maximum transmission and the in- 
put voltage is varied linearly using a triangular wave- 
form that varies between —Vt^ and This wave- 
form is generated by a fast function generator (Tektron- 
ics AFG3252) with two output channels whose frequen- 
cies and phases can be varied independently. The mod- 
ulated photons are then sent to SPCMs (Perkin Elmer 
SPCM-AQR-13), which are connected to the start and 
stop inputs of a TDC (Fast-Comtec TDC 7886S). Coin- 
cidence counts are binned into histograms and plotted as 
a function of the time difference between the detection of 
a Stokes and an anti-Stokes photon. 

In Fig. [3] (color online) we demonstrate the modulation 
of a biphoton wavefunction. The data are receded by bin- 
ning coincidence counts versus time into 1 ns bins. In Fig. 
|31Ja), the Stokes modulator is turned off and is biased at 
maximum transmission; and the anti-Stokes modulator is 
driven at 35 MHz. In agreement with theory, the bipho- 
ton wavefunction is not modulated. In FiglHJb), both 
modulators are modulated at 35 MHz with the same, but 
arbitrary, phase. The correlation function is now mod- 
ulated at twice the applied frequency, i.e. at 70 MHz. 
When driven by non-sinusoidal waveforms, the correla- 
tion function is modulated by the cross-correlation of the 
two modulating signals. In agreement with Eq. [3], we 
have verified that when the two modulating signals are 
square waves with the same frequency, the correlation 



function is modulated by a triangular function. 

We next demonstrate the Fourier measurement tech- 
nique. As shown in Fig. [2{a), we use a 35m long polariza- 
tion maintaining fiber to delay the anti-Stokes photon by 
175 ns. The modulators are driven synchronously at fre- 
quencies between MHz and 30 MHz. From the Fourier 
transform property that delay in the time domain corre- 
sponds to oscillation in the frequency domain, the plot 
of coincidence counts versus modulation frequency shows 
ripples with a frequency (5.7 MHz) equal to the inverse 
of the time delay. To simulate slow detectors we bin co- 
incidences into 1/xs bins, and since the temporal extent 
of the waveform is less than 1/is, the entire biphoton is 
contained within the first time bin. As a first example, 
we set the MOT optical depth to 15, the coupling laser 
Rabi frequency, Qc ~ 3.47r3i and the pump laser Rabi 
frequency, fip = 4.74r3i, where T31 = (27r) x 6 x lO*' is 
the spontaneous decay rate out of state |3). For these pa- 
rameters, the correlation function has an approximately 
Gaussian-like shape. In Fig. [D^a), we plot the number 
of coincidence counts in the first 1/is bin as a function 
of the modulation frequency. For this plot, the exper- 
imental rate of paired photons is 325 and for each 
point data is collected for 80s. In Fig. Hfb), we compare 
the biphoton wavefunction as directly measured with 1 
ns bins (black line) to the Fourier cosine transform of the 
trace of Fig. |4l(a) (gray line) vertically scaled for a one 
point fit. We find reasonable agreement between the two 
methods. We note that the sharp spike in the temporal 
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FIG. 5: (color online) Fourier Transform measurement tech- 
nique at an optical depth of 35. (a) Frequency domain data, 
(b) Fourier Cosine transform (grey), and temporal correlation 
(black). 



art polymer light modulators operating at 200 GHz [l9| 
and therefore modulating at 400 GHz, then the Fourier 
technique, will allow measurement down to about 2.5 ps. 
This is about eight times faster than the fastest reported 
superconducting detector [2^ . Looking further to the fu- 
ture, all-optical light modulators have been demonstrated 
at frequencies greater than 1 THz 2l|; thereby in prin- 
ciple allowing measurement of photons on femtosecond 
time scales. We note that an advantage of the Fourier 
technique as compared to either sum frequency correla- 
tion or Hong-Ou-Mandel interference is that the signal 
and idler photons may be correlated at distant detectors 
and do not need to be brought together at a summing 
crystal or at a beam splitter. 

The authors acknowledge helpful discussions with S. 
Sensarn and P. Kolchin. P. Kolchin also made substantial 
contributions to the apparatus used in this work. This 
work was supported by the U.S. Army Research OfRce, 
the U.S. Air Force OfRce of Scientic Research, and the 
Defense Advanced Research Projects Agency. 



trace near t ~ results from the DC component in the 
frequency domain trace. 

We next increase the optical depth to 35 and set 
ilc = 2.45F31 and rip = 3.2IF31. The biphoton wave- 
function is now characterized by a leading edge spike 
and an almost rectangular main body. FiglS^a) shows 
coincidence counts collected in the first 1/is bin plotted 
as a function of modulator frequency. Here, the experi- 
mental rate of paired counts is 80s^^ and for each point, 
data is collected for 200s. The cosine transform of this 
trace, multiplied by a vertical scaling factor, is plotted 
against time in Fig. Eljb) (gray line). For comparison, we 
show the correlation function where coincidences have 
been histogrammed into 1 ns bins. The shapes measured 
by the two techniques arc once again in reasonable agree- 
ment. 

In summary we have shown how biphotons may be 
modulated, and how this modulation may be used to 
measure the magnitude of the biphoton wavefunction. 
Though we have used long biphotons and low modula- 
tion frequencies, we believe that this technique should be 
extendable to short biphotons. A commercially available 
telecommunication modulator driven at a frequency of 
60 GHz and therefore modulating at 120 GHz [Eq. (4)] 
will allow measurement of biphotons with a minimum 
length of about 8 picoseconds. This is about a factor 
of five faster than state-of-the-art commercial SPCMs. 
If the comparison is made on the basis of state-of-the- 
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